NASA 

Technical Memorandum 107241 


Army Research Laboratory 
Technical Report ARL-TR-1 1 17 


A Computer Program for the Computation 
of Running Gear Temperatures Using 
Green’s Function 


S. Koshigoc, J.W. Murdock, and L.S. Akin 
Gearese arch Associates 
Banning, California 


D.P. Townsend 
Lewis Research Center 
Cleveland, Ohio 








■jg SB® 4 


CXI 



o> 

CXI 


CO 

o> 

o> 


Prepared for (die 

Seventh International Power Transmission and Gearing Conference 
sponsored by the American Society of Mechanical Engineers 
San Diego, California, October 6-9, 1996 



National Aeronautics and 
Space Administration 


U,S. ARMY 


m 

inrv. 


1 


RESEARCH LABORATORY 


T»istrebution statement a 

Approved for public release; 
Dialribuljon Unlimited 



A COMPUTER PROGRAM FOR THE COMPUTATION OF RUNNING GEAR 
TEMPERATURES USING GREEN S FUNCTION 


S. Konhkjoe/* J.W. Murdock/ and L. S. Akin*** 

Gearesearch Associates 
750 Indian Welts Rd. 

Banning f California 92220-5308 


D. P. Townsend*** 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland Ohio 44135 


ABSTRACT 

A new technique has been developed to study two dimensional heat 
transfer problems in gears. This technique consists of transforming the 
heat equation into a line integral equation with the use of Green’s theo- 
rem. Tlie equation is then expressed in terms of eigenfunctions that sat 
isfy the Helmholtz equation, and their corresponding eigenvalues for an 
arbitrarily shaped region of interest. The eigenfunction are obtained by 
solving an intergral equation. Once the eigenfunctions arc found, the 
temperature is expanded in terms of the eigenfunctions with unknown 
time dependent coefficients that can be solved by using RungoKuita 
methods. Ifiic time integration is extremely efficient. Therefore, any 
changes in the time dependent coefficients or source terms in the bound- 
ary' conditions do not impose a great computational burden on the user. 
The method is demonstrated by applying, it to a sample gear tooth. Tem- 
perature histories at representative surface local ons are given. 

INTRODUCTION 

The cooling of girars is an important problem that has been studied for 
a numlicr of years. An early model of oil cooling is given in Dc Winter 
and Block (1072). F.l-Bayourny cl al. (1980) expands on die model of 
De Winter and Block (1972) mainly by noting the importance of the 
Coriolis force on oil cooling of gears. LI- Bayou my et a! (1989) ;i)so 
develops a finite element model of the gear tooth. We now believe that 
the finite clement method Iras speed and accuracy limitations and have 
abandoned this approach in favor of the Green’s function method 
described herein, 

In die next section, a now dynamic, accurate and efficient solution 
method for two-dimensional heat transfer problems in gears is described. 
The solution method consists of transforming lire heat equation into a 
line integral equation with the use of Green’s functions with unknown 


time dependent coefficients. The transformed integral equation is used 
to obtain the dynamic equations for the time dependent coefficients for 
each eigenfunction. In order to obtain the eigenvalues and coirespond- 
mgeigenfu notion?;, the Helmholtz equation for the eigenfunctions is trans- 
formed into i\ line integral equation by the use of the two-dimensional 
free <;pacc Green’s function. The integral equation is discretized into a 
set of homogenous simultaneous equations. The discretized version of 
the eigenfunction can be obtained by solving a set of homogenous 
simultaneous equations. 'Flic obtained dynamic equations can tie inte- 
grated extremely efficiently. Therefore, any changes in t he boundary con 
riitinns do not impose a great computational burden. In the third section, 
the computational results and a discussion is presented. 

FORMULATION 

In this section, an accurate and efficient solution method for solving a 
imoe dependent two dimensional heat problems in gears is developed. 
The temperature field is expanded in terms of the eigenfunctions with 
unknown lime dependent dynamic coefficients. The dynamic equation 
for the time dependent coefficients of each eigenfunction is obtained by 
the use of an integral equation. 

Consider a gear tooth geometry shown in Fig 1 . The: transient heating 
of the gear tooth is cribed by 

^ = uV’T (i! 

dt 

whcic V 2 = J)'/i)x 2 + d 2 /Dy 2 and ot is the thermal diffusivily. The 
boundary condition for the left side and lop of the gear is given by 

K t= h M- T c) 
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where the equ alien represents the heat output to cooling oil. ^ is the 
vector coordinate on the boundary, K is the thermal conductivity, n is the 
outward pointing normal unit vector, h is the heat transfer coefficient 
andT c is the temperature of the cooling oil. The boundary condition for 
loaded side (meshing surface) of the gear (right side in Fig. 1) deter- 
mines the heat conduction into gear body and is given by 

Kg:.F 0 (t) (3, 

wheic F 0 is the beat flux into the gear and the remainder of the 
boundary ts described by the vanishing normal derivative ot die tem- 
perature. dT/Qn = 0. 


Dyna mical Equations fo r T ransien t Hmi Flow 

Rather than solve the phy sical cquarions directly, as in a finite dement 
method or finite difference method, we develop hero a Greco’s function 
method that reduces the two-dimensional problem to a onc-dimc ns in nal 
lino integral over the gear tooth boundary. The line integral equation 
yields the dynamical equations. Hi is procedure yields a computational 
advantage over the previous approaches where there is no reduction in 
the problem dimensionally and the time and spatial integrations are per- 
formed simultaneously. 



Figure 1. — Simulated gear tooth geometry a net* all the 
boundary conditions. 


We assume, the following eigenfunctions, v n . and eigenvalues, k n are 
knov-n and the eigenfunctions satisfy the Helmholtz equation, 

(V ? + k*)v»=0 (4) 

inside the gear tooth with vanishing normal derivative. dyjifo - 0. (The 
method for obtaining die eigerotiluo and eigenfunctions is given in the 
following section.) It can be shown that these eigenfunctions arc or 
thogonal. The eigenfunctioos will be normalised to unity hy 

JJviWy-i. 


The heat equation is transformed into a line integral equation by using 
Green's theorem. (Wyld, 1972) 



where ds is the surface element and dl is the line element. Substituting 
Fqs. (\) and (A) into the left hand side of Eq. (5) yields 



where the vectors a and c in die integration limits are shown in Fig. 1 
(a driven gearj. Note that a considerable simplification has been accom- 
plished for the line integral on the right hand side of Eq. (5). The only 
contribution comes from the boundary conditions described by the non- 
vanishing normal derivatives. Eqs- (2) and (3). Substitution of Eqs. (Z) 
and (3) and expansion of the temperature field in terms of the 
eigenfunctions as 


T(x.t) = ^n„(l)v n (x) 


(7' 


yield a dynamic equations for the time dependent coefficients.. r\ n . 


where 


t| n 4 CtkJ^ Tf n — + + Gf. 

L »■ 

A mn.=|. h (C)v«(I)Vn(5)dl 


(«) 


(9.0 


B 



(9b) 


C n =J>o(q ¥n (|)d. (9c) 

The line integrals sue easily calculated ai- the boundary conditions are 
changed. Therefore, (he calculation burden tor changing the boundary 
conditions is minimal. These coupled first order equations in Eq. (Ki are 
integrated efficiently with Runge-Kutta method 


Ei genvalues and eigenfun ction? 

Tlic method for obtaining the eigenvalues and eigenfunctions used to 
effect the simplification of the previous ‘icction is given below. These 
functions arc two dimensional analogs of the trigonometric func tions 
used in Fourier methods. 

The eigenvalues, fc and the corresponding eigenfunctions. \|/ n are 
obtained through the use of two dimensional free space Green s func- 
tion (Koshigot* andTubis 10H9). 
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( 10 ) 


* where H 0 f 1 5 is the Rankel function of the first bind arid the Green's func- 

tion satisfies the inhomogeneous He] mho Hz equation. 

; (v 2 J-ka)G(xix') ==-8(x- x') (11) 


where 5( x ) is the delta function. 

Again utilizing Green's theorem, one obtains 


( v ' 2 + k 2 ) V ;1 - V „ ( V' 2 4- k 2 )g) ds • 


■s 


■ 1 ' on 




dn' 


d!' 


H2) 


With die use of Eqs. (4) and (11) and vanishing normal derivative bound- 
ary condition for ihe eigenfunctions the integral equation is simplified 


f n (i)-Jv4-r)-^di'=o 


(I?) 


This integral equation yields d>; value of v n at any location inside the 
gear tooth when the values of \ji n on the gear tooth boundary arc known 
Now let x approach n point. £ on the gear tooth boundary, then Tiq. (13) 
becomes 


( II') 

P^v,(t)HPj v»(E')~^ — ;-<n' =0 < in 

Where is the contribution from die singularity in the integrand 
(Burton and Miller 197 1) and if. given by: 

= inside angle at the point q/2;t ( 15) 

and P designates the Cauchy pnncipal value integral The eigenfunction, 
y „ is discretized in Eq. ( 14) and yields a set of simultaneous equations. 
The eigenvalues, k n are determined by setting the determinant o( the 
simultaneous equations to zero. Once the 'dgciivalucft arc determined, 
the correspond), j eigenfunctions can be obtained through the 
simukiiiieuus equations. The formulation developed in tins section can 
now bo applied to the gear tooth geometry (shown in big 1) and the 
calculation result discussed. 

SAMPLE CALCULATION 

The technique developed in the previous section is applied to the gear 
tooth geometry shown in big. 1. Various coordinates are labeled in Pig. 1 


in order to specify the key features of the gear tooth geometry. The 
cooardinatc* of these points are given in table 1. The physical constants 
used for die calculation are: the gear thermal diffusivily, a = 0.452 fl 2 /hr; 
the thermal conductivity. K - 25 Btu/hr/ft/F; Die nil temperature, T c - 
200 F; the heat transfer coefficient, h -= 0.34 Btu/sec/ft 2 /F follow s from 
Dc Winter and Block < 1972) and El Baypoumy el al. ( 1 989). 

The heat flux, F 0 . along the boundary from the location £ = £ i.o c 
(shown in Fig. 1 ). is given as a function of the distance measured from 
die point £ (shown m Fig. 2). This is the heat generated in mesh fora 1 
inch wide gear, with a pitch radius of 6 inches, rotating a 10,000 rpm and 
transmitting 500 hp. We arc interested in gear steady-slate temperatures 
that take hundred*, of seconds to reach. Hence, the detailed tempenKure 
changes that occur as the gear goes in and out of mesh cannot be re- 
solved. Thus, both the \ <uing and the cooling have been averaged over a 
complete revolution cycle 


Table 1 . — Coordinate of 
Points in Figure 1 


Figure 1 

x, inches 

y, inches 

location 



a 

-0.197 

0 541 

b 

.073 

.795 

c 

.274 

.337 

d 

.3 IK 

.184 

c 

.371 

.184 

f 

.362 

.00(1 


The temperature calculation was performed u-'ing 21 cigenniodct. (This 
was found to provide adequate com ergcncc.l Thc re suits presented herein 
were generated on a 486 -PC with a total runniii): time, including the 
costly eigenfunction generation, oi less than one hour 
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Figure 2.— Time-averaged hont input to gear. 


.j 
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Figure 3. — Gear surface temperature at selected locations. 

llie eigenvalues and the corresponding eigenfunctions a le obtained 
through the use of discretized version of Eq. (14 ). The discretization is 
carried out with 40 line element:* and on each ek menu the second order 
approximation is used to represent the geometry as well as the 
eigenfunctions, 'the temperature calculations at three different locations 
on the gear tooth boundary. (0.0.0. 397) represented by the solid line, 
(-0.228,0.066) by the dashed line, and (0.174,0. 194) by the dotted line 
arc shown in Fig. 3 as functions of linn*. The curve labeled 1 1 is centered 
on rhe iop of the gear at location (0.0. 0.795). T2 on the cooled side of 
the gear in the region cooled by the oil at (-0.228.. 0.464). T3, the highest 
curve, is on the side healed during mesh, slightly above the pitch point at 
(0. 1 74, 0.592). Because both the heating and cooling functions Have been 
averaged over a complete gear revolution , these temperatures can best be 
interpreted as out-of-mesh temperatures. 

At the present Lime the calculation method is hard wired for a single, 
but representative, f problem. Furthermore, we have not yet included fthe 
triangular portion of the gear extending to the axis of rotation. Generali- 
zation of the method is planned now that it* application to a specific 
problem has been demonstrated. Appendix A "Them) a! Analysis of Spur 
Gears" provides the geometric formulation input. 

CONCLUSION 

A new technique has been developed to study two -dimension a? heat- 
ing nf gears, The computational advantage of this technique over 


previous approaches using the linite element method or the finite dif - 
ference method results from two features: First, the problem is reduced 
from two dimensions to one. Second, the lime and spatial integrations 
arc separated. Therefore, when compared will) other methods, this new 
technique can provide substantial improvements (one order of magDi 
tude) in computational speed. However, the benefit of the dimensional 
reduction is manifested not only in the computation speed but also in 
the ease of problem set-up since one is required to deal with the bound- 
ary not the entire two-dimensional gear geometry. The other benefit of 
this technique, based on the r eparation of time and spatial integration, 
is accuracy. This technique lakes the full advantage of the spectral 
method that has exponential solution convergence. This should be com 
pared to finite element or difference methods where only algebraic con- 
vergence is possible. 

Since the new technique is extremely efficient, any changes in the 
time dependent coefficients or source terms in the boundary conditions 
do not impose a groat computational burden on the user. This result is 
very important when performing accurate Scoring Analysis in gears. 
Tliis uI lows l c bulk (or blank) temperature to be accurately known. 
Furthermore, the gear out-of-mesh temperature is not a constant along 
the tooth profile at steady state running conditions, as is often assumed 
by gear engineers. The method is also more adaptable for use in small, 
lubricated, concentrated contacts, such as gears, since high resolution 
can be obtained without using large numbers of elements. 

Currently, wc are planning to extend this technique to a multiblock 
application that includes the remainder of the gear sector and that funner 
optimizes the computation accuracy and speed. 
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APPENDIX A 

THERMAL ANALYSIS GEOMETRY FOR SPUR GEARS 

by 

Lee S. Akin * and Dennis P. Townsend ** 


to the gear temperature analysis problem. The above is accomplished by 
lining special integration techniques newly developed to accommodate 
ibe special conditions found in lubricated concentrated sliding contacts, 
such as found in high performance aircraft gear drives. This appendix 
describes the intricate gear geometry needed for this analysis, mostly 
available in the literature but certainly not all in one place This paper 
hopes to fill that need. 


ABSTRACT 

The gear geometry needed to perf onn a complete thermal analysis for 
a gear set is very' complex and has, to the authors knowledge, never been 
published before a*, one set of equations needed for the analysis. The 
thermal analysis of gears is a very important subject in that it can be used 
to determine the scoring failure criteria which includes the blank (hulk) 
temperature and the flash temperature as used in the Blok scoring 
formula 


I T, = T, + Tv, 

The geometry for thxs analysis must include the involute geometry, 
the load variance when (he teeth mesh is tn a sinije or double tooth 
contact 7X)nc, tooth load sharing due to varying deflection, long and 
short addendums. Up and/or root relief. Also included are the affeci of 
heat partitioning due to varying too th contact sliding velocities over con- 
stantly changing Hertzxan contact band widths, gear sc I speed and its 
affect on tlic lubrication regime (film thickness versus surface rough- 
ness). 

'["his geometry analysis will be used as the input parameters to com- 
plete the analytical computation of running gear temperatures using 
Green's function. 


INTRODUCTION 

This is one more in a series of papers by the authors on iheir continued 
study of the ari of the prediction of the onset of scoring or scuffing in 
high performance gear drive mechanisms and our attempt to make it 
more scientific. Those studies have examined a scrief: of disciplines from 
Jn interdisciplinary lubrication theory (ref. I), m a study of the cfiect of 
windage on the lubricant flow into high speed gear teeth (Ref. 2) and a 
model for lubricant flow in between gear teeth (Refs. 3. 4, 6, 7. id 
and 11, 13, 16). 

Jn addition these studies have evaluated die analytical and experimen- 
tal spur gear temperature effects on opeialing variables (Ref. 5). gear 
lubrication and cooling studies {Refs. 1 to !), arid 12), investigations of 
affect of transient (time variant) thermal and lubricant boundary layers 
(Refs. 14 and 15). 

Most of the abm * and all of the thermal work was done using finite 
element methods which produce large matrixes causing slow computer 
solutions to provide satisfactory accuracy. This geometry analysis is 
developed for, "a computer prop am for the compute linn of running gear 
temperatures using Greens function” provides a new and unique solurion 


FORMULATION OF GEOMETRICAL EQUATIONS 

The equations for the involute etuve in rectanyulaj coordinates is shown 
in equation set 2 (*cc big. 1 ). 


y - Involute curve 



2.x-- (cosf) -r 9sin0j R or x R sinO see nomcncaliuic 
y - (sin 9 ! 0cos9jR h or y=RcosO 

x : *y 2 =0 + e 2 ) l ' 2 -h 

Kh 


* Gcarescarch Associates. ** NASA-Lowii Research Center. 
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The poiar/vector/il angle 0, in Frig. 1 may he calculated from: 


/ i ? n l/2 1 

> = (x-+y J -R b ) 


(R 2 -Rg) 


I/? 


-tan 


(r 2 -r l) 


V2 


Rb 

~ tan<fr-4 ,=5 inv6 


where: R = --&R 0 - ~AR b = Rcm< 
so that She j>res<.ure angle 4> in Jig. 1 is: 


(R 2 *Rb)' /2 ,r k 

4. G -- tan 1 ;= cos'* 1 -—and the roll angle is f. = 0 r 4> = 

R b R 

tan$. The important radii ol curvature equations. some at critical loca- 
tions* axe shown below in equation .‘el 5, ami sec Figs. 1 and 2: 

5. p = fcsinp = eR b and for die mating gear p H1 ~ R (1| sini|) = E m R bm 
I “he lowest point of contact for the mesh is p c calculated from 

p c =Csin«-(R|+Rf;) 

Pcm ~ (Rom + R b: 


\l'-' 


and 


on 


* bin I 


Rr == 


- 2 

^Csin»-VR5 + Rbl ■ ,? b I 

where C sin $ - L a the line -of action, and 


il/2 


^Csin^-^R^ -#-K blll ) 


,// 2 


0 G 



Figure 2. — Gear mesh line ol action- geometry. 


?ri) 

is. P-— The distance of roll/slide S along the involute curve 
may be calculated by integrating over the roll angle e, from (vogle. 
Ref. 17), sec Fig. 1: 


is the lowest point on the mating ge Jr. The radii of curvature at the low- 
est and highest points of single tooth contact may be calculated from: 

pi = (Ro* K b ) 1 " Ri=(pi -Rb) 1 '* 

/ ) \t ( 2. 

Ph’Csm^R^-R^,) +p b 

R h -(Pb *- Ku)" 

2n 

where the base circle pitch is Pb and the circular pitch 


6. 




and since, as can be noted from Fig. 1, p- R b f, = R b Lin<J> tuule^— • 
2 2 R * 

tfius S|_n - — over any arbitrary portion o / the profile with sub- 

2R b 

2 2 
c Pn “Pc 

scripts 1 and 2 and -o ” over Die who]*- profile from the? fir«i 

point of mesh contact on the profile at $ c to the last poml at outside 
diameter S fl . For example: fiom R c to Kj to R (at pitch point) to R h to 
R 0 al the outside radius. The arc length along the base circle £cnn he 
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calculated f ror»« ^were subscripts 1 and 2 arc at arbitrary locations on the 
tooth profile) a;; shown. 

7. 2 — ) ~ R b ■*" t(inOj)so that X j _ 2 = 

^ t ‘( R~ ~ R^" j = ^ *" Pi and thus the ratio between the distance S and 

the distance £ Incomes k == a quantity useful in thermal calcu- 

2R b 

labors. The transition time, as a function of radii of curvature and its 

robing velocity between any two points along the tooth, can be expressed 
in equation set 8 as: 


8 . 



dr rad 

--R h - R b u> — 
dt sec 


where: 


dc (j) do 

*T' - 6) , a constant end since p = fR k then dr. = <adt and — dl = — , so 
dt D f n 


that 


fpl 


dp. , p; 


Pi 


0) / , to(f^-q) 

(/j — jf| ) — — — i-i and At. , 

e >- c ' i (fc-.,) 


(to-t,) 


P2-P1 


In ” 


R b .ji p, 

'ideal dimoi 
ure 2 and described! below in equation 9: 


'$ the fame it takes for the gear to rotate from Cj 


to £3. 'Hie critical dimensions aJoiig the line* of action are shown in fig- 


9. Z - V R h - R b j + V^o2 - r ^2 -CsinA is the length of the "line of 

contact b- f as a subset of the ‘Mine of action a - g 11 

Now’ we can calculate the *viddi of the Hertzian band of mutual con- 
tact at the inesh point from equation 10 (Timoshenko Ref.l 8; 



Figure 3.— Gear tooth profile showing normal load W^. 


Iliercforc the rolling velocity anywhere along die profile ear. be cal- 


culated from: V 


np. , v - r ; 

— t-j — ft/sec and V ir> - 77—7:7 ft/sec 


114.59 


114.59 


for die mat- 


ing gear so that the sliding velocity is calculated from v s 


H P, + n mPtfn 
1 14.59 


fi/scc where n = speed of gear and n in - the speed of the mating gear 
Anotiier value needed to calculate die coefficient of friction is the total 


np + n m p :ry . 

velocity from: v t -jjj- — ft/scc. 

We can now calculate the coefficient of friction as: 


r ( 3.17x 10 6 W. scriji 


where V - tooth face width. 


Pep ” viscosity in cemipoisc and W t = tangential tooth load in lbs. Thus 
we fiwd ourrelves in a position to calculate the instantaneous heat flux 
qfp)j us a function of the radius of curvature at the instantaneous posi- 
tion along Uk line of contact per equation Set 13. 


13. 


t i(p.)= 


W, scc(6)t(np, -n m p, ro ) 
170 330 


BTU/nain 


( 16 W n (K 1 tK 2 )pi1 >2 ^ 

1 F(Pi + Pd ) 


where: 



and 


- l “ y 3 _ V 1 - v 1 - I 

^ : ~ anrf K 1 + K 2 ~ n £ I and since: + p 2 s L a == C riinjj 

and W n r: W f /cosri = VV^e.c^, see Fig. 3, and E, = = : E 


B = 3.19 


\FcosPC4in.pE- 


l/.i 


Tlic rolling velocity for the gear V, and h's mating gear V 2 arc calcu- 
lated from equation set II: 

11. v l “ “TrrJr 1 and v 2 2! ‘^7^' ft/scc. so that the sliding 

30£) 3oU 

velocity V f - V| - V 2 ■= (njpj -n 2 p; ) it/*ec when;- Vj = V' & = 

Vm. 


Ihis equation can also be written in a form more useful using the 
pinion speed only: 


<i(p,)=- 


W, sec($)n p f 


170 330 


"(Ppi “"'gpgi) BlU/min, where 



the 


gear ratio At times it is more convenient to calculate the hear flux from 
tile radius vector at the instantaneous contact points where: 


qR, 




foT R pc _ 

M K f“ W 


2 V 

S l |u r i fir 

v /2 


-R 


S pV 




ps J 


BTU/min 


CLOSURE 

This appendix develops the geometry analysis for \hc input for 
the computer solution of line thermal analysis of spur gears using 
Green s function to solve for gear blank surface temperatures. This 
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geometry analysis provides the gear tooth geometry input, the 
gear tooth sliding distance parameters, the rolling and sliding 
velocity inputs and the equations for the frictional heating 
developed during the gear tooth meshing, as a function of the 
location on the gear tooth. 

Using these inputs the program can then determine the transient 
and steady state temperatures of the gear teeth. 
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NOMENCLATURE FOR APPENDIX A 


R 


O 


B 

C 


width of the band of mutual contact at the mesh point 
center distance for mating pair 


R,&R h 


outside radius of the gear and end of the involute 
curve and R orn for mating gear 

radius vector to lowest and highest points of single 
tooth contact 


E 


Young's modulus of elasticity 


V 


velocity along the curve (involute) 


Y 

f 

M* 

N 

n 

P 

Ph 

<i R i 


face width 

coefficient of friction 
gear ratio N^/N p 


V 5 the sliding velocity at an instantaneous point Vj - v, 

in the mesh 

V t the total velocity at an instantaneous, point in the mesh 


number of Iccth in gear and N (J) number of teeth in 
the mating gear 

pinion speed, rpnu n m for mating gear and n p 
for pinion 

diametral pitch of the teeth 


x & y 


rolling velocity of the gear and its mating gear 

tangential and normal (perpendicular) load, 
respectively 

Cartesian coordinates of the involute curve from its 
origin at the base circle 


circular pitch of gear and its mating gear Z 

heat flux calculated using the instantaneous radius 

vector At j 2 


length of the line of contact as a subset on the line 
of action 

the time it takes to rotate from Ej to r 2 




V* ‘v 


r s . ar p. 




1PC1 


heat flux due to sliding friction at instantaneous point 
of contact using the radius of curvature as a 
parameter 

radius vector to the pitch point and R m lor the mate 
at the contact point 

base radius of the involute curve (its origin) and 
for mating gear 

radius vector to lowest point of contact from center 
of gear and R crtl for mating gear 

base radius of gear and pinion respectively 

instantaneous ntdius og gear and pinion respectively 

outside radius of gear and pinion respectively 

lowest point of contact for gear and pinion 
respectively 


AN 

C 

6 

P 

Pc 


Pi & Ph 

* 


virtual number of teeth expansion or reduction for 
long and short addendums 

roll angle on tooth 

involute polar angle = inv6 = tan$ - <b - £ 0 (rad' 

oil viscosity, in cp 

radius of curvature from the base circle and p m for 
the mating gear at contact point 

radius of curvature tit lowest or initial point of 
contact from base circle and p crp for mating gear at 
contact point 

radius of curvature at lowest and highest points of 
single tooth contact 

pressure angle of mesh 
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